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ABSTRACT 

This work presents the main ultraviolet (UV) and far-infrared (FIR) proper- 
ties of two samples of nearby galaxies selected from the GALEX (A = 2315A, 
hereafter NUV) and IRAS (A = 60/mi) surveys respectively. They are built in or- 
der to get detection at both wavelengths for most of the galaxies. Star formation 
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rate (SFR) estimators based on the UV and FIR emissions are compared. Sys- 
tematic differences are found between the SFR estimators for individual galaxies 
based on the NUV fluxes corrected for dust attenuation and on the total IR lumi- 
nosity. A combined estimator based on NUV and IR luminosities seems to be the 
best proxy over the whole range of values of SFR. Although both samples present 
similar average values of the birthrate parameter b, their star-formation-related 
properties are substantially different: NUV-selected galaxies tend to show larger 
values of b for lower masses, SFRs and dust attenuations, supporting previous 
scenarios for the star formation history (SFH). Conversely, about 20% of the 
FIR-selected galaxies show high values of 6, SFR and NUV attenuation. These 
galaxies, most of them being LIRGs and ULIRGs, break down the downsizing 
picture for the SFH, however their relative contribution per unit volume is small 
in the local Universe. Finally, the cosmic SFR density of the local Universe is 
estimated in a consistent way from the NUV and IR luminosities. 

Subject headings: surveys: GALEX — ultraviolet: galaxies — infrared: galaxies 



1. Introduction 

What is the best way to measure the SFR of galaxies on large scales and at different 
redshifts? The possibility of estimating the SFR of a galaxy directly from the luminosity at 
a single wavelength would be a major advantage for anyone wanting to compute the SFR per 
unit volume at a given redshift. This quantity could be derived directly from the luminosity 
function (LF) at this wavelength and at this redshift. Under these conditions, large area 
surveys at single wavelengths might suffice. 

The recent SFR of a galaxy is often measured from the light emitted by the young stars: 
given their short lifetimes their luminosity is directly proportional to the rate at which they 
are currently forming. The UV and FIR luminosities of star forming galaxies are both closely 
related to the recent star formation: most of the UV photons are originally emitted by stars 
younger than ~ 10 8 yr, but many of these photons are re-processed by the dust present in 
galaxies and re-emitted at FIR wavelengths. Strictly speaking, neither of these fluxes can 
be used alone to estimate the SFR independently of the other one (e.g. Buat & Xu 1996; 
Hirashita et al. 2003; Iglesias-Paramo et al. 2004). Because of the previous lack of data 
at both wavelengths, attempts have been made using only the rest frame UV (Lilly et al. 
1996; Madau et al. 1996; Steidel et al. 1999; and more recently Schiminovich et al. 2005 
with GALEX data), or just FIR data (Rowan- Robinson et al. 1997; Chary & Elbaz 2001), 
but only a few authors have compared both (Flores et al. 1999, Cardiel et al. 2003). The 
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SFR estimator based on the UV luminosity suffers from attenuation by dust and it has to 
be corrected in order to properly trace the SFR: for UV-selected samples of galaxies the 
attenuation can reach more than 1 mag (e.g. Iglesias-Paramo et al. 2004; Buat et al. 2005). 
On the other hand the FIR emission is not free of problems because the dust can also be 
heated by old stars, and can be a non-negligible correction for many star forming galaxies 
(Lonsdale & Helou 1987; Sauvage & Thuan 1992). Neither of these two indicators taken 
alone is an accurate estimator of the SFR except perhaps for starburst galaxies where (a) 
the dust attenuation was found to follow a tight relation with the slope of the spectrum at 
UV wavelengths (Meurer et al. 1999), thus allowing one to estimate the dust attenuation 
with only information on UV fluxes, and (b) the contribution to the dust emission coming 
from old stars can be neglected (Sauvage & Thuan 1992). In the most general case, the 
best estimator of the SFR should contain combined information of the luminosities at both 
wavelength ranges (Hirashita et al. 2003). The UV and FIR fluxes are thus complementary 
for tracing star formation and it is well known that the FIR/UV ratio is a proper indicator 
of the dust attenuation (Buat et al. 1999; Witt & Gordon 2000; Panuzzo et al. 2003). 
Although other indicators of the recent SFR of galaxies have been extensively used in the 
literature, a detailed discussion of their quality as SFR tracers will not be discussed in this 
paper. 

The GALEX mission (Martin et al. 2005a) is imaging the high-Galactic latitude sky at 
two UV wavelengths (A = 1530A, FUV; A = 2315A, NUV) and is providing the astronomical 
community with unprecedented data (both in quantity and quality). The UV data combined 
with existing FIR datasets (from the IRAS, ISO or Spitzer missions) now allow us to carry 
out detailed studies of the UV and FIR properties of galaxies, with special emphasis on the 
derivation of the dust attenuation and star formation activity in star forming galaxies. 

With this purpose in mind we selected two samples of galaxies: one from the GALEX- 
All Imaging Sky (AIS) and the other from the IRAS PSCz (redshifts, infrared and optical 
photometry, and additional information for 18,351 IRAS sources, mostly selected from the 
Point Source Catalog) and FSC (Faint Source Catalog) for which UV and FIR fluxes are 
available. With these datasets in hand we undertaken a study of the properties related to 
their emission at these wavelengths. Both samples were extracted from the same region of 
the sky (~ 600 degrees 2 , constrained by the status of the GALEX survey when this work 
was initiated). From the GALEX catalog we built a complete sample of galaxies down to 
AB NUV = 16 mag 1 and cross-correlated it with the IRAS database (FSC), allowing non- 
detections at 60/xm (fluxes lower than 0.2 Jy). The FIR-selected sample was built from the 



1 AB magnitudes are denned as: AB V = — 2.5 log /„ — 48.6, where f v is the monochromatic flux density 
expressed in erg s _1 cm~ 2 Hz -1 . 
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IRAS catalog (PSCz) with a limiting flux at 60/xm of 0.6 Jy as the only constraint. The 
resulting list of PSCz sources was cross-correlated with the GALEX database, allowing again 
non-detections in the NUV. Both the NUV and the 60/zm limits used to build the samples 
were chosen to allow for a very small number of non-detections and to sample the galaxy 
population over a large range of values of the dust attenuation. Besides an analysis of the 
SFR, for these NUV and FIR-selected samples of galaxies (chosen with well-defined selection 
criteria) can also be used to place important constraints on models designed to predict the 
statistical properties of galaxy populations. The attenuation of star light by interstellar dust, 
and its emission in the far infrared are usually computed very crudely in models of galactic 
evolution. Dust attenuation is usually deduced in such approaches from other quantities 
such as the mass gas and the metallicity (e.g. Guiderdonni & Rocca-Volmerange 1987; 
Devriendt & Guiderdoni 2000; Balland et al. 2003). The properties of large samples of 
galaxies observed both in the UV and FIR with clear selection criteria, such as the ones 
presented in this paper, provide an important statistical constraint for the calibration of the 
treatment of dust in such models. 

The first paper in this series (Buat et al. 2005) based on these samples was mainly 
devoted to the dust attenuation properties. In the present work we discuss various aspects 
relating to the NUV and 60/xm emission of our sample galaxies, including systematic differ- 
ences in the 60/im and NUV luminosities and we will focus on their star formation related 
properties. This paper is organized as follows: the samples are presented in section 2. The 
relation between NUV and 60/zm luminosities is discussed in section 3. Section 4 is devoted 
to the derivation of the SFR and to a comparison of various estimators as well as to a discus- 
sion on the star formation activity related properties of the samples. The derivation of the 
local cosmic SFR density by different methods is discussed in section 5. The main conclu- 
sions are presented in section 6. Throughout this paper we adopt the following cosmological 
parameter set: (h,Vl ,\ ) = (0.7,0.3,0.7), where h = H /100 (kms~ 1 Mpc _1 ). 



2. Observational dataset 

2.1. NUV-selected sample 

In order to build the NUV-selected sample, we used 833 fields of the GALEX All sky 
Imaging Survey (AIS) each having an exposure time equal to or larger than 50s. We used 
only the central 0.5 deg radius circles in each field to ensure a uniform image quality: the 
resulting sky coverage is 615 deg 2 . Within this area we selected all the galaxies of the 
GALEX AIS survey with AB NUV < 16 mag. This bright limit was chosen in order to ensure 
IRAS detections of all the galaxies with attenuation larger than ~ 0.3 mag (for a limit of 
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0.2 Jy at 60 /iin from the IRAS FSC using the calibration of Buat et al. 2005) and highly 
significant upper limits for the less attenuated galaxies. 

On the one hand, the moderate angular resolution of GALEX (FWHM ~ 6 arcsecs) 
does not allow for a secure discrimination of stars from small galaxies. On the other hand 
the GALEX pipeline can induce some shredding of larger galaxies. This results in multiple 
detections which (cumulatively) correspond to a single galaxy because of the misidentification 
of star forming regions as if they were individual galaxies. The main consequence of this 
is the underestimation of the fluxes of large galaxies. Corollary catalogs were thus required 
in order to perform a reliable selection of galaxies. Our starting point was the catalog 
of NUV detections produced by the GALEX pipeline 2 , which made use of the Sextractor 
code (Bertin & Arnouts 1996). We made the assumption that all the galaxies brighter than 
ABnuv =16 mag, even if they are shredded, should contain at least one Sextractor detection 
brighter than AB NUV = 18 mag. Then we associated an object from databases of well 
known stars and galaxies (SIMBAD, 2MASS, LEDA) with each of the Sextractor detections 
brighter than AB^uv — 18 mag. The problem of shredding was mostly resolved by using 
LEDA. As this database contains the optical diameters of the galaxies, NUV detections 
corresponding to shredded galaxies can be associated with their parent galaxies provided 
they are located within the aperture defined by the optical diameters and the position angle 
listed in LEDA. For the detections of galaxies not shredded we used SIMBAD and 2MASS 
in order to classify them as stars or galaxies. Finally, objects not associated with any known 
source were classified as "dubious" . In order to test the quality of our selection method, we 
cross-correlated our final catalog with the SDSS DR1, which spatially overlap one fifth of 
our sample. The spectral and photometric information of the SDSS together with its higher 
angular resolution made possible an optimal classification of all the objects detected in the 
field: all the objects present in both GALEX and SDSS catalogs were found to be properly 
classified. Dubious objects were found to be noise detections or ghosts generated mainly near 
the edges of the GALEX frames. Thus we ended up with a catalog of bona fide galaxies or 
fragments of galaxies (i.e. belonging to the large galaxies shredded by Sextractor) brighter 
than ABnuv = 18 mag. The next step was to obtain aperture photometry of these objects 
in order to select out only the galaxies brighter than AB NUV = 16 mag. 

Photometry of our sample galaxies was performed in the GALEX NUV and FUV bands. 
Since the selection criterion for our sample galaxies was imposed in the NUV band, we also 
took the NUV images as our reference for the total photometry. We performed aperture 
photometry using a set of elliptical apertures, the total photometry corresponding to the 



2 version 0.2.0, September 2003, with the correction to the NUV and FUV magnitudes reported in Buat 
et al. 2005 
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aperture where convergence of the growth curve is achieved. Once we determined the to- 
tal NUV flux, the photometry in the FUV band was obtained by performing background 
subtracted aperture photometry within the same elliptical aperture where convergence of 
the NUV growth curve was achieved. This way, the NUV— FUV colors are consistent. Some 
galaxies were contaminated by the flux of nearby bright stars or galaxies. The contaminating 
sources were then masked in the NUV and FUV frames in order to obtain proper NUV and 
FUV fluxes for our galaxies. Table 1 shows the typical uncertainties of the NUV and FUV 
magnitudes. NUV and FUV magnitudes were corrected for Galactic extinction using the 
Schlegel et al. (1998) dust map and the Cardelli et al. (1989) extinction curve. In the end, 
a total of 95 non-stellar objects brighter than AB NUV = 16 mag were found. One of them, 
classified as a QSO, was excluded of the sample. 

FIR fluxes at QOfim were obtained from the IRAS Faint Source Catalog (FSC, Moshir 
et al. 2000) for 68 of our 95 galaxies. We discarded all these sources for which the cirrus 
parameter (as listed by the FSC) is larger than 2 because it can result in uncertain fluxes. 
The Scan Processing and Integration Facility (SCANPI) was used to obtain the FIR fluxes for 
the remaining 27 objects. Three of these galaxies (UGC 11866, UGCA 438 and UGC 12613) 
were not detected at 60/zm. We adopted a conservative upper limit of 0.2Jy at 60/im (as 
given in the FSC) for these galaxies. Four galaxies of the sample were not covered by the 
IRAS survey. 

2.2. FIR-selected sample 

The FIR-selected sample was extracted from the IRAS PSCz (Saunders et al. 2000) 
over the 509 deg 2 in common with the GALEX AIS fields having exposure times larger 
than 90 sec. In order to keep only good quality FIR data we discarded those galaxies for 
which the probability of a correct optical identification of the FIR-selected galaxies was lower 
than 50%, as listed in the PSCz. As for the NUV-selected sample, galaxies for which the 
cirrus parameter (as listed in the PSCz) was larger than 2 were discarded. A total of 163 
galaxies were selected; all but two of them (Q00443+1038 and Q23367-0448) have published 
radial velocities. As galaxies were selected from the PSCz, the imposed limiting flux at 
60/im was 0.6 Jy. This limit, combined with the one estimated for the GALEX AIS at NUV 
(ABnjjv ~ 20.5 mag, Morrissey et al. 2005), results in detections at NUV for all the galaxies 
with dust attenuation as large as ~ 4.4 mag (Buat et al. 2005). Indeed, only two galaxies 
(Q00443+1038 and Q00544+0347) were not detected in the NUV frames and a total of 23 
were not detected in the FUV frames. 

The NUV and FUV photometry of the FIR-selected galaxies was performed using the 
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same technique as for the NUV-selected galaxies. 

2.3. Completeness and bias of the samples 

Before drawing conclusions about the properties of the samples we have to check on just 
how representative they are. Because of the reduced statistics of the samples, if the sampled 
volumes are not large enough it could be that some luminosity ranges are oversampled or 
undersampled with respect to reference samples defined over larger volumes of the local 
Universe. 

We check the representative nature of our samples by building LFs and comparing them 
to the standard ones at z — 0, constructed from larger samples of galaxies (NUV LF of Wyder 
et al. (2005) and 60//m LF of Takeuchi et al. (2003)). We calculate both NUV and 60/im 
LFs of our sample by the Lynden-Bell method (Lynden-Bell 1971), implemented to obtain 
the normalization using the formulation of Takeuchi et al. (2000). The calculation of the 
uncertainty is based on a bootstrap resampling method (Takeuchi et al. (2000)). We note 
that the Lynden-Bell method is robust against density inhomogeneities, and hence we can 
trust the LF so determined (see Takeuchi et al. 2000 for details). The results are shown in 
Figure 1. The error bars correspond to la uncertainties. The agreement between the LFs of 
our samples and the corresponding LFs from larger samples is very good, so we are confident 
that in spite of their small size our samples are representative of flux-limited NUV and FIR 
samples in the local Universe. 

We also compare the volumes from which the samples were extracted. Figure 2a shows 
the redshift distributions for both samples. The median values are 0.013 and 0.027 for the 
NUV and FIR-selected samples respectively. At a first glance, this means that the FIR 
selection samples a volume 8 times larger than the corresponding NUV selection. However, 
we must recall that the redshift distribution of a flux limited sample is strongly dependent 
on the shape of the LF, and as shown in Figure 1, the NUV and 60fim LFs are very differ- 
ent. In Figure 2b we show the theoretical redshift distributions for NUV and FIR-selected 
samples with the same limiting fluxes as our two samples 3 . As can be seen, the theoretical 
median values of the redshift for both samples are consistent with the ones obtained from 
the observational data. A limiting magnitude of AB NUV = 18 mag is required to obtain 
similar median values of the redshift distributions of both samples, as we show in Figure 2c. 
And in this case most of the galaxies detected in NUV will not have any counterpart in 
FIR. This behavior of the redshift distribution can be understood intuitively. Indeed, the 



3 Dctails on the calculation are given in Appendix A. 
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flux density selection procedure omits intrinsically low-luminosity objects from the sample, 
whereas bright objects are hardly affected by the flux selection. To the degree that the LF 
is well reconstructed from the flux/magnitude-limited samples, these samples can be said to 
be representative, with respect to the luminosity and/or flux density, and it is indeed the 
case for the present work. 

3. Relation between L 60 and L NUV 

The physical link between the FIR and UV luminosities of galaxies is rather complex. 
On the one hand, both are related to the light of young stars, so one expects a correspondence 
between the two. On the other hand, the FIR emission is due to the absorption of the UV 
light thereby leading to an anti-correlation. Since our samples were built to avoid upper 
limits - i.e. most of the galaxies selected in NUV (or at 60 /im) are also detected at 60 /xm 
(or in NUV) - we are able to discuss statistically the intrinsic relation between both the two 
wavelengths and outline the specifics of NUV versus FIR selection effects. 

In Figure 3 we plot L^yy versus L^q for both samples 4 . The two samples exhibit 
very different behaviors: the NUV-selected galaxies show a good correlation between both 
luminosities, with a correlation coefficient (in log units) of r ~ 0.8. On the contrary the 
dispersion is very large for FIR-selected galaxies and the correlation coefficient is low: r ~ 
0.3. 

NUV-selected galaxies appear intrinsically less luminous at 60/iin than FIR-selected 
ones. This is also true for the sum of both luminosities, L tot = L NUV + L 60 , which is 
supposed to be a crude estimate of the bolometric luminosity of galaxies related to recent 
star formation (e.g. Martin et al. 2005b). Although the luminosity distribution within each 
sample is the combined effect of the LFs and selection criteria, this result is confirmed by 
other studies and from the comparison of the 60fim and NUV LFs themselves (e.g. Martin 
et al. 2005b; Buat et al. 2005). Both distributions flatten at higher luminosities, reflecting 
a general increase of the dust attenuation already pointed out in the literature by several 
authors (Wang & Heckman 1996; Buat et al. 1999; Sullivan et al. 2001; Vijh et al. 2003). 

One could argue that the difference in luminosity between the two samples is a con- 
sequence of bias in the sampling. We show in Figure 3 the lines corresponding to our 
lower (upper) limit of the NUV attenuation above (below) which the NUV (FIR) selected 



4 Throughout the paper the NUV and 60 fim luminosities Lmuv an d L§ will be calculated as vL l 
expressed in solar units. The adopted value for the bolometric solar luminosity is L & = 3.83 x 10 33 erg s _1 
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sample is complete. Thus, the fact that only very few low-luminosity and low-attenuation 
FIR-selected galaxies are detected must be taken as real. Low-luminosity, high-attenuation 
galaxies should have been detected by our FIR survey if they were present. For the same 
reason, very luminous galaxies should have been detected in the NUV survey if they existed. 

The good correlation found between L NUV and L 60 for the NUV-selected galaxies has 
to be related to their rather low dust attenuation: in these galaxies both L NUV and L 60 
represent a significant part of the total luminosity of the galaxies. This result holds for 
intrinsically faint galaxies (L tot < 2 x 1O 9 L ). The very loose correlation found for FIR- 
selected galaxies may also be explained by the effect of the dust attenuation. With a mean 
dust attenuation larger than 2 mag, the NUV luminosity becomes a poor tracer of their total 
luminosity whereas the 60/im luminosity is not very different from the bolometric emission of 
the young stars. Some fluctuations in the percentage of NUV photons escaping the galaxies 
can induce large variations in the NUV observed luminosity on an absolute scale without 
any strong physical difference on the scale of the total luminosity of the galaxies. 

We make a final comment on the so-called "UV luminous galaxies" (UVLGs, defined 
as those with L FUV > 2 x 10 10 L Q in Heckman et al. 2005). We found 3 ULVGs in our 
NUV-selected sample and a total of 8 (including the previous 3) in the FIR-selected sample. 
All but one of these galaxies are more luminous at 60fim than in the NUV (in fact most 
of them are LIRGs), and their attenuation is typically larger than 1 mag. This means that 
these galaxies are not only UV luminous but also very luminous from a bolometric point of 
view. 



4. Selection effects on observational quantities and physical properties of 

galaxies 

The main aim of this section is to show the effect of the selection criteria of samples 
of galaxies on observational and physical quantities. We will now show that the selection 
criteria of a sample of galaxies play an important role in defining the nature of the galaxies 
selected and thus, in their averaged properties. Accordingly we warn against the unqualified 
comparison of results obtained from samples of galaxies selected on the basis of different 
criteria. 

In order to reduce the uncertainties associated with the FIR and NUV fluxes we impose 
further constraints on our galaxy sample: 

• Ellipticals, SOs as well as AGNs (Seyferts and QSOs) were excluded since the origin 
of their 60/zm and NUV fluxes is clearly not associated to recent star formation. The 
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necessdary classification information is available for most of the NUV-selected galaxies; 
but this turned out not to be the case for the FIR-selected galaxies, so contamination of 
the sample by ellipticals and/or AGN among these galaxies cannot be totally excluded. 
Galaxies with extraneous radio sources (from NVSS and/or FIRST) within the IRAS 
beam were also excluded since part of the FIR flux of these galaxies could be due to 
contaminating background objects. 

• Multiple galaxies, not resolved by the IRAS beam but clearly resolved into various 
components in the GALEX frames were excluded, since a one-to-one 60/im-NUV as- 
sociation is not possible for them. 



After applying these criteria we ended up with 59 and 116 galaxies from the original NUV 
and FIR-selected samples, respectively. Hereafter we will use these restricted subsamples for 
our subsequent analysis of the star formation related properties, although we will keep the 
terminology FIR and NUV-selected samples to refer to the restricted subsamples. Given 
that all the galaxies were extracted from the same region of the sky, some of them belong to 
both subsamples. Their basic properties are listed in Table 2: (1) Identifier of the galaxy; (2) 
"Y" ("N") indicates whether the galaxies is included (or not) in the NUV-selected sample; 

(3) "Y" ("N") indicates whether the galaxies is included (or not) in the FIR-selected sample; 

(4) R.A.(J2000 equinox) of the source ; (5) Declination (J2000 equinox); (6) Radial velocity 
in km s _1 , obtained from NED or LEDA; (7) Distance to the source in Mpc, corrected for 
the Local Group Infall to Virgo and H = 70 km s -1 Mpc -1 ; (8) Morphological type, from 
NED or LEDA; (9) IRAS identifier: "F" for FSC; "Q" , "O" or "R" for PSCz; "SCANPI" 
for absence in both catalogs. 

Table 3 gives some useful photometric data for the galaxies in the restricted subsamples: 
(1) Optical identifier; (2) NUV magnitude corrected for Galactic extinction; (3) FUV mag- 
nitude corrected for Galactic extinction; (4) Flux density at QOfim in Jy; (5) Flux density at 
lOOyum in Jy; (6) H magnitude from 2MASS Extended Source Catalog. For galaxies with 
no detection by 2MASS we adopt the limiting value of H = 13.9 mag (3 mJy) as given by 
Jarrett et al. (2000); (7) NUV attenuation in mag, derived as in Buat et al. (2005); (8) 
FUV attenuation derived as indicated in Buat et al. (2005). For some galaxies Eqs. 1 and 
2 gave negative values of the NUV and FUV attenuations which is, of course, unphysical. 
In fact, this is an artifact of the polynomial fitting used to derive a functional form for the 
attenuation in Buat et al. (2005). Throughout this paper they will be considered as zero. 
Given that the FIR fluxes were extracted from different catalogs (PSCz for the FIR-selected 
sample and FSC/SCANPI for the NUV-selected sample), for those galaxies belonging to 
both samples we list the FIR entries corresponding to the PSCz catalog. For those galaxies 
not present in the FSC and not detected by SCANPI at 60/xm, we list f 60 < 0.2Jy, which 
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is the nominal limiting flux of the FSC. No estimate of an upper limit at 100 jum is given 
for these galaxies. For galaxies with no detection in 2MASS we adopt the limiting value of 
H = 13.9 mag (3 mJy) as published by Jarrett et al. (2000). 

In Table 4 we list some star-formation properties which will be used in the forthcoming 
discussion: (1) Identifier of the galaxy; (2) SFR NUV from Eq. 3; (3) SFR FUV from Eq. 4; 
(4) SFR dust from Eq. 5; (5) SFR tot (NUV) from Eq. 6; (6) SFR tot (FUV) from Eq. 6 but 
modified by using SFRp UV instead of SFR^uy] (7) (SFR) averaged over the galaxy's 
lifetime, estimated as indicated in Appendix B. 

4.1. SFR derivations 

This section is devoted to a detailed comparison of the recent SFR as seen in the FIR 
and NUV-selected samples. Although other estimators of the recent SFR can be found in 
the literature (see Kennicutt 1998 for an interesting review on several methods to derive 
the SFR), we focus on only two of them, those using the NUV and FIR fluxes. Our aim 
is to compare commonly used recipes to derive SFR from the UV and FIR luminosity of 
the galaxies. Therefore we will make very classical calculations, as described below. For 
consistency we re-derive the calibrations in a homogeneous way, adapted to the GALEX 
bands: the formulae are found to be very similar to those of Kennicutt (1998). 

The underlying physical justification for deriving the SFR of a galaxy from the UV 
luminosity is the following: most of the UV photons emerging from a galaxy originate in 
the atmospheres of stars younger than ~ 10 8 yr. Thus, the SFR is proportional to the UV 
luminosity emitted by the young stars under the assumption that the SRF is approximately 
constant over this timescale. This is reasonable given that Salim et al. (2005) and Burgarella 
et al. (2005) found that the intensity of the youngest burst in large samples of nearby galaxies 
contributes typically less than 5% to the total. However, the presence of dust absorbs a part 
of the UV light escaping from galaxies and breaks down the proportionality between the 
SFR and the observed UV luminosity. 

As star-forming galaxies may present a large variety of relative geometries between 
stars and dust, the scattering of the stellar photons through the interstellar medium may 
introduce a fraction of them in the line of sight before they escape the galaxy. Thus, the 
effect of the dust differs from a pure extinction but is a complex combination of absorption 
and scattering. Following Gordon et al. (1997) we will use the term 'dust attenuation' for 
this global process at work in galaxies. The most commonly accepted method to estimate 
the dust attenuation at UV wavelengths is to use the ratio of FIR-to-UV fluxes (Buat & Xu 
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1996; Meurer et al. 1999; Gordon et al. 2000). Several analytical expressions are already- 
available in the literature for different UV wavelengths (Panuzzo et al. 2003; Kong et al. 
2004; Buat et al. 2005). All these expressions are fairly consistent except at high values of 
the dust attenuation, where some dispersion appears (e.g. Meurer et al. 1999, Kong et al. 
2004, Buat et al. 2005 at A ~ 1500A.). In this work we use the prescription of Buat et al. 
(2005) to obtain the corrected NUV and FUV luminosities: 

A NUV = -0.0495x 3 + 0.4718a; 2 + 0.8998a; + 0.2269 (1) 
where x = log Ljr/ L^uv and 

Apuv = -0.0333y 3 + 0.3522?/ 2 + 1.1960?/ + 0.4967 (2) 

where y = \ogL IR /L FUV . 

Once the observed NUV and FUV luminosities have been corrected for dust attenuation, 
the SFRs can be derived using the following expressions 5 : 

log SFR NUV (M & yr- 1 ) = log 

J^NUVjCorr (L )-9.33 (3) 

\ogSFR FUV {M & yr~ l )=\og 

-^FUV,corr (L e ) ~ 9.51 (4) 

In Figure 4 we show the ratio of SFR NUV / SFR FUV as a function of L tot (= L NUV +L eo ), 
which traces the bolometric luminosity related to recent star formation and has the advantage 
of being a purely observational quantity. As this figure shows, both quantities are equivalent 
with a dispersion of about 20%. Since our sample is NUV-selected, hereafter we will use 
NUV as our reference wavelength for star formation related properties. 

The luminosity at IR wavelengths provides a different avenue to the derivation of the 
SFR. Dust absorbs photons at UV wavelengths and re-emits most of them at IR wavelengths 
(8 — 1000/im). Under the hypothesis that all the UV photons are absorbed by dust, the IR 
luminosity would be a direct tracer of the SFR of a galaxy. One source of uncertainty is 
the difficulty in estimating the total IR luminosity from the FIR flux at only one or two 
wavelengths. In this paper we use the prescription of Dale et al. (2001) and derive L IR by 
using / 60 and / 100 . For the galaxies for which only / 60 is available we use the mean value of 
iW/ioo estimated using the galaxies detected at both wavelengths. If we assume the same 
scenario as for Eq. 3, the SFR can be expressed as: 

log SFR dust (M Q yr- 1 ) = log L IR (L Q ) - 9.75 (5) 



5 This formula has been derived from Starburst99 (Leitherer et al. 1999) and assuming solar metallicity, 
and a Salpeter IMF from 0.1 to 100 M . 
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However, Eq. 5 is a good approximation only for the most extreme starbursts, since many 
of the FIR-selected galaxies are, in fact, detected at UV wavelengths. A further limitation 
of this method concerns the fraction of the total IR luminosity heated by old stars (the 
cirrus component, hereafter represented by 77), which should be removed before applying 
Eq. 5. This quantity is known to depend on the morphological type of galaxies (Sauvage & 
Thuan 1992), but a precise estimate for individual galaxies is subject to large uncertainties 
(Bell 2003). 

The SFRs estimated from these methods are often compared in the literature for indi- 
vidual objects or for large samples of galaxies. In order to see whether they are consistent 
with each other we show here a comparison of the two using the galaxies of our two sam- 
ples. Figure 5 shows the ratio SFR NUV / SFR dust as a function of L tot ; each sample shows 
a different behavior. For the NUV-selected sample (blue, filled circles), SFR NUV is always 
larger than SFRdust (and the ratio can be as high as 3) but the discrepancy is lowered as 
L tot (and A^uv) increase. This result is expected since we have seen in Section 3 that low 
luminous galaxies are brighter in the NUV than at 60/im. This affirms that SFRdust cannot 
give a proper estimation of the SFR for these galaxies. 

The FIR-selected galaxies extend the trend found for the NUV-selected sample to higher 
luminosities. For values of L to t > 3 x 10 10 (and for higher values of the dust attenuation), 
where no NUV-selected galaxies are present, SFRdust systematically exceeds SFR^uv by 
a factor of ~ 2. One reason that could play a role in this inconsistency between the two 
estimators is that the dust attenuation is not properly estimated for very dusty galaxies. In 
any case it does not make sense to use the corrected UV luminosity to measure the SFR for 
these IR bright galaxies. In fact, Charmandaris et al. (2004) have reported decoupled IR 
and UV emissions for some dusty galaxies, which could be at the basis of the discrepancy 
found between SFR NUV and SFR dust found in this work for galaxies with large attenuation. 

The conclusion of this analysis seems to be that SFR^uv is a good tracer of the SFR 
for low values of the attenuation, and in the opposite extreme SFRdust must be used for very 
heavily attenuated galaxies. There is no obvious way to delimit these two different regimes, 
or to chose which and which of the two indicators should be used in the intermediate cases. 
And so we warn users about any undiscriminated comparison of SFR NUV and SFR dust for 
samples of galaxies selected with different criteria. 

An alternative tracer of the SFR containing information from NUV and IR luminosities 
has already been discussed in the literature (Hirashita et al. 2003; Iglesias-Paramo et al. 
2004, Bell 2003): 
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SFR tot = SFR° NUV + (1 - 77) x SFR dust (6) 

where 77 accounts for the IR cirrus emission and SFR° NUV is obtained following Eq. 3 but 
using L^uv,obs (that is the observed NUV luminosity) instead of L^uv,corr- This estimator 
has the advantages of being free of the model dependence of the attenuation correction, and 
it contains information of the observed NUV and the IR luminosities. 

One limitation of this estimator, 77, is the adopted value of the IR cirrus contribution. 
Hirashita et al. (2003) and Iglesias-Paramo et al. (2004) reported a value of 77 ~ 0.4 for 
normal disk galaxies. Accurate values of 77 for individual galaxies are not easily obtained 
and instead, averaged values are often used. However, this parameter is strongly dependent 
on the sample of galaxies under study and cannot be easily generalized. Whereas an average 
value of 77 ~ 0.4 seems to apply for normal disk galaxies, a value of 77 ~ seems to better 
represent the properties of starbursts (Hirashita et al. 2003). Bell (2003) also proposed 
a cirrus correction for a compilation of galaxies from the literature with FUV, optical, IR 
and radio luminosities. He found 77 ~ 0.32 ± 0.16 for galaxies with L IR < 10 11 L and 
77 ~ 0.09 ± 0.05 for galaxies with L IR > 10 11 L . For our NUV-selected sample (similar 
to the normal star forming galaxies of Hirashita et al.) a value of 77 ~ 0.2 gives similar 
values for SFRnuv and SFR to t(NUV). Although our NUV-selected sample must contain 
galaxies more active than that of Hirashita et al. (since their selection is based on optical 
fluxes rather than on UV fluxes), this result gives an idea of the uncertainties related to the 
determination of 77. For practical issues, throughout this paper we will use the value of 77 
of Bell (2003) - not far from that of Hirashita et al. (2003) - when computing SFR tot , but 
keeping in mind that the uncertainties reported by this author are of the order of 50%. 

Another limitation of SFR tot is that it depends on the wavelength at which we measure 
the UV flux. In order to illustrate this point we show in Figure 6 the ratio of SFR tot (NUV) / SFRt ot (FUV) 
as a function of SFR tot (NUV) for both samples. As can be seen, for the NUV-selected 
galaxies SFR tot (NUV) is systematically larger than SFR tot (FUV) by about 20%. This 
discordance for the NUV-selected galaxies is due to the fact that the UV attenuation is not 
grey: A NUV < A FUV for most galaxies (see Buat et al. (2005) and Table 3), and since we 
showed in Figure 4 that SFR NUV m SFR FUV , it is obvious that SFR Q NUV > SFR FUV . On 
the contrary, for the brightest FIR-selected galaxies the agreement between SFR tot (NUV) 
and SFR to t{FUV) is good since for these galaxies SFR to t is dominated by SFRdust- We 
conclude that SFR tot is stable to within 20% for whatever UV wavelength at which we 
measure the UV flux. 

We compare now SFR^Qf^NUV) to the classical estimators S F Rnuv and S F Rdusti in 
order to set their domain of applicability. Figure 7a shows the comparison between SFRnjjv 
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and SFRtot- At low values of the SFR both quantities are almost identical for the NUV- 
selected galaxies. This is expected since for these galaxies both A^uv an d Lj R are almost 
negligible and SFR NUV ps SFR° nuv . As the SFR grows, we note an increase of SFR NUV 
with respect to SFR tot , but always within ~ 15%. This increase could be due to the choice 
for the cirrus correction and/or to the fact that A^uv does not exactly corresponds to the 
dust emission (since factors other than absorption do play a role in the attenuation, like for 
example the relative geometry between stars and dust). Finally, the NUV-selected galaxies 
with the largest values of SFR show a decrease of SFR NUV with respect to SFR tot . We 
stress that these galaxies have L IR > 10 11 L and so their cirrus correction is different 
from for the rest. All in all we find that for the NUV-selected galaxies, basically those 
with SFR NUV < 15 M & yr -1 , SFR NUV and SFR tot are equivalent to within ~ 15%. The 
FIR-selected galaxies show a different behavior. Whereas those with Lj R < 10 11 L Q show 
an ~ 15% excess of SFR NUV with respect to SFR tot , similar to the NUV-selected galaxies, 
for those with Lj R > 10 11 L , SFR^uv is well below SFR tot . This is easily understood 
as a consequence of the already mentioned discrepancy between SFRdust and SFR^uv for 
galaxies dominated by their IR emission. 

In Figure 7b we compare SFRdust and SFR tot . The NUV-selected galaxies follow a 
very dispersed trend with SFRdust I 'SFRtot increasing with SFR to t- This behavior is due to 
the fact that SFR dust lacks the UV contribution which is dominant in these galaxies. The 
FIR-selected galaxies obey two different trends: for galaxies with SFR tot < 15 M Q yr~ x 
the ratio log S F R dust / SFR tot ~ 0, although with a dispersion of ~ 0.2 dex. This large 
dispersion is due to the contribution of the NUV luminosity to SFR tot , which is important 
for the less attenuated galaxies. On the contrary, at large values of SFR to t the average value 
of log S F Rdust I ' S F Rtot ~ 0.04 dex with a very small dispersion. This is a consequence of the 
fact that most of these galaxies have Lj R > 10 11 and are dominated by their IR emission, 
so the difference between SFRd us t and SFR to t corresponds basically to the cirrus correction 
applied to SFR to t, which is minimal. 

Bell (2003) proposed a calibration of the SFR similar to the one described in Eq. 6 but 
using FUV as the reference UV wavelength. His method is based on the relation he found 
between L IR /L FUV and L IR (L IR / L FUV ~ \J Lj R /10 9 ) for a compilation of galaxies from 
the literature with FUV, optical, IR and radio luminosities. One can see in Figure 8 that 
our NUV-selected galaxies follow well the Bell's relation whereas it is not the case for the 
FIR-selected sample. The galaxy sample used by Bell is therefore closer to a UV selection 
than to an IR one. Again we emphasize the importance of the selection biases in deriving 
SFRs. 

The overall conclusion emerging from this study is that SFR tot seems to be a proper esti- 
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mator of the SFR of galaxies whatever their dust content is, since it avoids the main problems 
of the clasical estimators SFR^uv an d SFRd us t and is consistent with them within their 
respective domains of applicability to within ~ 15%. We again warn against indiscriminate 
comparisons of the SFR of galaxies estimated from these classical estimators since the results 
could be strongly affected by selection biases as we have illustrated in this section. The com- 
bined uncertainty of SFR tot due to the choice of the UV wavelength at which we measure 
the UV flux and to the cirrus contribution to the IR luminosity is < 55%. Throughout this 
paper, we adopt SFR tot (NUV) as our proxy to trace the recent SFR. 



4.2. Star formation history 

The determination of the SFR of a galaxy gives information about the total number 
of young stars that are being formed. But this does not necessarily mean that the light 
coming from this galaxy is dominated by these young stars given that most galaxies are 
composed of a mixture of various stellar populations of different ages. This parameter is of 
major importance in understanding the SFH of the Universe. Recent results based on large 
amounts of SDSS data suggest that the higher the mass of a galaxy, the earlier its stars were 
formed (Heavens et al. 2004), thus supporting the so-called "downsizing" explanation for the 
SFH of galaxies already proposed by several authors (e.g. Cowie et al. 1996; Brinchmann 
& Ellis 2000; Boselli et al. 2001). We devote this section to the study of the SFH of the 
galaxies in our samples. 

A quantitative estimation of the SFH of a galaxy requires information relating the 
relative contribution of young and old stars. The birthrate parameter (hereafter b) has been 
proposed as a quantitative estimator of the recent SFH of a galaxy (Scalo 1986). It is defined 
as the ratio between the current and the past-averaged SFR: 



Since b depends on the overall SFH of the galaxy, an accurate estimation from obser- 
vational quantities is complex and involves several approximations. A detailed derivation 
of b following the prescriptions of Boselli et al. (2001) can be found in Appendix B. As 
explained in Section 4.1, the NUV luminosity is sensitive to the SFR over a timescale of 
~ 10 s yr, and thus b is not sensitive to shorter-timescale variations in the SFH. However, 
this is not a serious problem since Burgarella et al. (2005) have shown that less than 20% 
of the galaxies in either sample have bursts younger than 10 8 yr. 

Figure 9 shows the distributions of b for both samples of galaxies. The median values 
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of both distributions are similar: 0.50 and 0.58 for the NUV and FIR-selected galaxies 
respectively. In Figure 10a we show the relation between SFR tot and b for both samples. 
In the range of overlap between the two samples (approximately 0.5 < log SFR tot < 1.5) 
the values of b are consistent, but beyond this region two different trends are seen: the 
NUV-selected galaxies show no trend of b with the SFR, whereas the FIR-selected galaxies 
show an increase of b for high SFR. This bimodal behavior of b is also seen in Figure 10b, 
where b is plotted as a function of the attenuation. Again, for the NUV-selected galaxies b at 
lower values of the attenuation, although this trend is very dispersed. The opposite holds for 
the FIR-selected galaxies, with galaxies with high b being the most attenuated. Thus, the 
picture emerging from this study is that galaxies dominated by young stellar populations fall 
into two categories: those showing low SFRs and low attenuation, which naturally appear 
in UV surveys, and those with high SFRs and large attenuation, mainly detected in FIR 
surveys. 

4.3. The link between the H luminosity and star formation properties of 

galaxies 

The baryonic mass of galaxies is a key parameter in understanding their evolution. It 
has been proposed as the parameter which governs the SFH, rather than the morphological 
type, for example (Boselli et al. 2001). In addition, it is often used to derive some properties 
like the dust attenuation in semi-empirical models of formation and evolution of galaxies. 
For this reason we devote this section to a discussion of the effects of the sample selection 
on the relation between the mass and the star formation related properties of galaxies. As 
explained in the previous section, we will use the if-band luminosity as a tracer of the galaxy 
mass. 

First we show in Figure 11 SFR tot as a function of the if -band luminosity. Both sam- 
ples show a positive relation between these two quantities, which means that more massive 
galaxies are also currently forming more young stars. This result is expected since we are 
comparing two extensive quantities. However, whereas the relation followed by the NUV- 
selected galaxies shows a small dispersion, the FIR-selected galaxies exhibit a more dispersed 
relation, especially at the most massive end. At high galaxian masses the range in SFR spans 
almost two orders of magnitude, which is not seen in the NUV-selected sample. 

In Figure 12 we show the dust attenuation as a function of the if-band luminosity. 
Two different trends are seen in this plot. The NUV-selected galaxies show a fairly good 
correlation between the two quantities, with the dispersion increasing towards high if-band 
luminosities. On the contrary, the FIR-selected galaxies span an interval of almost 5 mag 
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in dust attenuation and no correlation at all is shown with the galaxian mass. While the 
trend followed by the NUV-selected galaxies could be interpreted as a result of the mass - 
metallicity relation reported for samples of spiral and irregular galaxies (Garnet & Shields 
1987; Zaritsky 1993) in the sense that more metallic galaxies contain more dust, there is no 
simple explanation for the lack of any trend shown by the FIR-selected galaxies. 

Finally, we show in Figure 13 the b parameter as a function of the iJ-band luminosity. 
The NUV-selected galaxies follow the classical trend that low-luminosity galaxies have larger 
values of b (e.g. Boselli et al. 2001). Some of the FIR-selected galaxies also follow this trend, 
although about 20% of them that show large masses and large values of b. As shown in 
Figure 12, these galaxies are among the most attenuated of the FIR-selected sample. Overall, 
our galaxies are shifted towards higher values of b with respect to the sample of galaxies of 
Boselli et al. (2001). These authors adopt a slightly different IMF than we do (M up = 80 M Q 
against 100 M ) and different evolutionary synthesis codes. Nevertheless, the large shift in 
b found between the samples can probably not be explained by these differences alone. The 
correction for dust attenuation could also partially explain the shift in b, since Boselli et al. 
assume average values of 0.20 mag for Sds and later types, and 0.60 mag for types earlier than 
Sd. For our NUV and FIR-selected samples, the values of the dust attenuation estimated 
from the FIR/UV flux ratio of each of the individual galaxies show higher averaged values 
for the two categories of morphologies than those of Boselli et al., which would imply higher 
SFRs. However, this effect is diluted by the fact that for many objects in their sample, 
Boselli et al. estimate the SFR as the mean value of SFR Ua and SFR UV . A further factor 
that could be responsible for the shift in b is the different selection effect of the sample: the 
sample of Boselli et al. (2001) is drawn from the nearby clusters Virgo, Cancer, Coma and 
A1367 and from the Coma-Al367 supercluster. Although not a unique selection criterion 
was applied, this sample can be defined as an optically selected sample of galaxies with a 
normal Hi content. Thus, in their sample there is a non-negligible fraction of Sa-Sab, bulge- 
dominated galaxies, which tend to lower the average value of b (see their Figure 2). Since 
our selections are based on NUV and FIR fluxes, we argue that we are surely avoiding these 
kind of objects. Anyway, one important point is that the NUV-selected galaxies follow the 
same relation between mass and b as the optically selected ones (disregarding the absolute 
calibration of both quantities) and that a fraction of the FIR-selected galaxies do not follow 
this trend. 

We have seen that the relation of the star-formation-related properties with the mass of 
galaxies strongly depends on the selection procedure of the sample: whereas for NUV-selected 
galaxies low-luminosity galaxies are also low mass, show low attenuation and have high values 
of b. A selection based on the FIR fluxes yields a different result: a population having high 
attenuation, high mass and strong star-formation activity appears. This population is absent 
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in the NUV-selected sample. Since these galaxies present very high values of the attenuation 
(most of them are LIRGs and/or ULIRGs), their UV (and optical) fluxes are strongly dimmed 
and for this reason they are often excluded from flux limited surveys. However, even if these 
galaxies show such extreme properties, they do not put into question the downsizing picture 
for the SFH of galaxies since their contribution to the local cosmic SFR density is very low 
(see Takeuchi et al. 2005). 



5. The local cosmic SFR density from different estimators 

We saw in Section 4.1 that a proper estimation of the SFR is not possible with informa- 
tion restricted only to either NUV or FIR fluxes. However, big surveys usually provide large 
amounts of data only at single wavelengths, thus an estimation of the density of SFR over 
cosmological volumes has to be carried out under these constraints. The accuracy of the 
cosmic density of SFR has already been discussed by Hirashita et al. (2003) using FOCA 
UV data. In this section we make a similar analysis using the new GALEX data. 

The usual way to estimate the average SFR density is to construct the monochromatic 
LF and then to weight the corresponding contribution to the SFR at a given luminosity with 
the probability of finding a galaxy with this luminosity. 



Px = k x J L x 0(L x )dL x (8) 

where p x is the SFR density estimated from the flux at a given A, k x is the conversion factor 
between SFR X and L x , and 4>(L X ) is the LF. 

A simple calculation using Eqs. 3, 8 and the NUV LF of Wyder et al. (2005) yields a 
value of Pnuv = 0.0091^^1 M © yr -1 M P C ~ 3 for the local cosmic SFR density. Unfortunately, 
there is not a simple relation linking the observed NUV luminosity and the attenuation (see 
Figure 3 in Buat et al. 2005), so we adopt the median attenuation of our NUV-selected 
sample, which is A^uv — 0.78 mag. After correcting for this median attenuation we obtain 
PNUV,corr = 0.018^ 3 M yr" 1 Mpc~ 3 . 

Takeuchi et al. (2005) report a value of the local cosmic star formation density derived 
from L IR of pm = 0.013+q;qq5 M yr _1 Mpc~ 3 . However, we recall that this quantity does not 
account for the fraction of UV escaping photons or for the cirrus IR heating. For our FIR- 
selected sample, the median contribution of the NUV escaping luminosity to SFR to t is 17% 
which leads to p' IR = O.OlG+Q^Jg yr" 1 Mpc~ 3 . After correcting for the cirrus contribution 
assuming rj = 0.32 we obtain pi^ corr = O.OIIIq;^ M yr _1 Mpc~ 3 which is well below 
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PNUv,corr, although still within the 1-a uncertainty. Although we have previously accepted 
that a cirrus correction of rj = 0.32 is valid for galaxies with Li R < 10 11 L Q and r] = 0.09 
for brighter galaxies, we argue that adopting just rj = 0.32 for all the galaxies is not a bad 
approximation for this calculation since it can be seen in Takeuchi et al. (2005) that the 
contribution of Li R (f>(Li R ) to the total J Li R (f)(Li R )dLi R of galaxies with Li R > 10 11 L is 
very low and can hardly affect our calculations. 

Proceeding in an analogous way as in Section 4.1, we can add both contributions to 
get the total cosmic SFR, and we get p to t = Pnuv + (1 — v) x Pm ( = 0.009 + 0.009) = 
0.018 M Q yr" 1 Mpc -3 , which is in good agreement with pNuv,corr, which means that correct- 
ing pnuv with a median attenuation is a good approximation. We also stress that the p to t 
is almost equally shared between the NUV and the IR contribution. 

The discrepancy found with pi R alone may be due to the average corrections applied. 
In order to obtain a better agreement between pNuv,corr, Ptot and pm tC orr two points should 
be studied in more detail: 

• A detailed knowledge of the cirrus component is required since this contribution is 
probably multivalued for a given value of L IR . Although the morphological type seems 
to drive this parameter (Sauvage & Thuan 1992), it could also be dependent on b 
since this parameter also measures the relative weight of the young and old stellar 
populations. A more detailed study of a large samples of galaxies, covering a wide 
range of galaxian properties, could shed light on its fractional contribution to the total 
cosmic L IR density. 

• The bivariate LF 4>(L NUV) L IR ) appears to be the best way to estimate the fraction 
of UV photons escaping from the galaxy, required to correct pi R . It is also required 
since for large values of L^uV) the attenuation can take on multiple values, and thus 
an average value, as the one used in this work, might be not the most appropriate. 

Under these conditions, the cosmic SFR densities would be expressed as: 



PNUV,corr = / / KNW X L NUV X IQ^nuv^uv ,L IR )/2.5 x 0( Lw> L IR )dL NUV dL IR (9) 




and 



PlR,corr 



II 



k ir x [1 - r)(L IR , L nuv )]xL ir x(1+(k nuv / ki R ) x (Law/ L IR ))x(j)(L NUV , L IR )dL NUV dL IR 

(10) 
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or if we use the approximation of Eq. 6 




Ptot — I I (k>ir x [1 — v(Lir, Lnuv)\ x Ljr + knuv x Lnuv) x 4>(Lnuv, LiRjdL^uydLiR 

(11) 



6. Conclusions 



We have performed a detailed study of the star formation properties of two samples 
of galaxies selected on the basis of their NUV and FIR fluxes, which were found to be 
representative of the nearby Universe when compared to samples drawn from larger volumes. 
The main conclusions of this work are: 

1. L NUV and L eo are tightly correlated for the NUV-selected galaxies. The opposite holds 
for FIR-selected galaxies, which span a large range of L 60 for a given value of L NUV 
and show larger values of attenuation. Intrinsically bright galaxies are more luminous 
at FIR than at NUV wavelengths, including the UV luminous galaxies (UVLGs), and 
they show moderate to high attenuation. 

2. The SFR deduced from the NUV fluxes, corrected for the dust attenuation (SFR NUV ), 
are not found to be consistent with those calculated using the total dust emission 
(SFRdust)- Whereas SFR^uv is larger than SFR^ust for galaxies with low attenuation 
(Anuv ^ 1 mag) the inverse is found for bright, but highly extinguished galaxies, 
mostly selected in IR: SFR^jjy is likely to underestimate the actual SFR in these 
galaxies by a factor ~ 2. A combined estimator based on UV and IR luminosities with 
a cirrus correction depending on the IR luminosity seems to be the best proxy over 
the whole range of values of SFR. As a practical recipe we found that SFR tot and 
SFR NUV yield similar results for SFR tot < 15 M yr _1 , whereas SFR tot and SFR dust 
are almost equivalent for SFR tot > 15 M yr _1 . 

3. NUV-selected galaxies follow a trend whereby low-mass galaxies show lower SFRs, 
low attenuation and higher values of b, indicating the existence of a dominant young 
stellar population. On the contrary, about 20% of the FIR-selected sample shows high 
attenuation, high SFRs and also large values of b, most of them being LIRGs and/or 
ULIRGs. In spite of their discordant properties, these galaxies are not sufficiently 
abundant in the local Universe to question the downsizing picture for the SFH seen at 
z = from optical surveys. 



4. The cosmic SFR densities of the local Universe, estimated from the NUV and IR 
luminosities, are consistent to within l-cr, although the difference between the two 
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values is large, when average corrections for the attenuation, UV escaping photons 
and IR cirrus component are applied. The sum of the individual contributions is quite 
consistent with the value obtained from the NUV luminosities corrected for average 
attenuation. A better knowledge of the cirrus contribution to Ljr and of the bivariate 
LF is required in order to better understand the large differences found between the 
monochromatic estimators of the local SFR density. 
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A. Relation of the Mean Redshift of Galaxies and the LF: Representativity of 

a Survey Depth 



We show the strong dependence of the mean redshift on the shape of the LF. This 
means that the mean distance of galaxies does not represent the depth of a survey, but 
rather reflects an intrinsic property of the sample. 

Since we treat a local sample of galaxies, we first approximate their distance by the 
classical Hubble's law: 



cz 



r 



(Al) 



where c is the velocity of light, and r represents the distance. We define M to be surface 
density of galaxies on the sky, and denote the LF as an explicit function of the characteristic 
luminosity, as 0(L/L*). Then, Af is written as (Peebles 1993): 



Af 




L(r) 



dL{r) 2 



r dr . 



(A2) 



-23- 



Using the detected flux density S, this can be expressed as 
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We observe the distribution function of the sources with a fixed flux density S as 
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The mean redshift of a flux-limited survey (limiting flux density S), (z) > s, is then 
defined as 6 



(z)>s = 



oo poo j2 



_ JS JO 



d 2 M 

dzdS' 



zdzdS' 



oo poo r2 



n 

Js Jo 



d 2 M 
dzdS 1 



(A6) 



dzdS' 



The numerator is obtained as 

y°° i d 2 M 



dzdS 1 



zdzdS' 



2a 2 \H 



1 



2a 2 \H 



4>{x)x 2 dx 



(j){x)x 2 dx 



S'~ 3 dS' 



s 

s- 2 



(A7) 



Here, x = az 2 S', a luminosity normalized by the characteristic luminosity of the LF, and 
expressed in terms of flux density S'. Since the luminosity must be positive, the lower bound 
on the integration with respect to x is 0, and upper bound is large, effectively taken to be 
+oo. Similarly the denominator becomes 
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6 Note that this is a different quantity defined by Equation (5.134) of Peebles (1993). 
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Both the numerator and the denominator of this part is the moment of the LF. This means 
that this is a function of its shape. 

Combining Equations (A6), (A7), and (A8), we have 



Let us carefully examine Equation (A9). First, the dependence of the mean redshift on 
the limiting flux density S is a power of —1/2. Also, it has the same order of dependence 
on L*. In contrast, the integral part of Equation (A9) has an important meaning. Since 
this part depends on the second-order moment, the tail of the LF affects the value very 
strongly. As we mentioned, we integrate over the (normalized) luminosity up to a certain 
very large value, and this part is a ratio between the moments of order 3/2 and 2. Hence 
the contribution from a large value of x controls the value. Consequently, the mean redshift 
(z) > s is very sensitive to the LF shape. 

This aspect is clearly seen in the comparison of the expected redshift distributions in 
NUV and 60/im calculated from the LFs and the limiting flux density or magnitude, because 
the shapes of the LFs at these wavelengths are very different (Buat & Burgarella 1998; 
Takeuchi et al. 2005). We show the comparison in Figures 2b and 2c. In these figures we fix 
the limiting flux density at 60yum to be 0.6 Jy, while we change the limiting magnitude from 
AB NUV = 16.0 mag (the actual value we adopt in this work) to 18.0 mag. The medians for 
the two wavelengths are very different when we adopt a limiting value of AB^uv = 16.0 mag. 
Only when we adopt a limiting value of AB^uv = 18.0 mag, which corresponds to a very 
sensitive survey, the median values approach for the NUV and FIR-selected samples. 



Boselli et al. (2001) give a detailed recipe to estimate (SFR), based only on observable 
quantities and/or adopted values for the parameters (see Gavazzi et al. 1996). Here we 
follow their prescriptions, using the if -band luminosity to estimate the past averaged SFR 




(A9) 



B. Estimating the birthrate parameter 
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and the same parameters as Boselli et al. (2001): 

_ gFg x t x (1 - R) 

L H x (M tot /L H ) x DM cont { ' 

where SFR in this work is averaged over 10 8 yr, t is the age of the disk (assumed to be 
equal to 13 Gyr), R is the fraction of gas re-injected by stars through stellar winds into the 
interstellar medium during their lifetime (taken to be equal to 0.3 for a Salpeter IMF), Lh 
is the if -band luminosity estimated as \og Lh = 11.36 — 0.5 x H + 2 x log-D (in solar units) 
where D is the distance to the source (in Mpc), M tot is the dynamical mass at the B-band 
25 mag arcsec -2 isophotal radius, M tot /L H is taken to be constant and equal to 4.6 (in solar 
units) and DM cont is the dark matter contribution to the M tot / Lh ratio at the optical radius, 
assumed to be equal to 0.5. 

REFERENCES 

Balland, C, Devriendt, J. E. G., & Silk, J. 2003, MNRAS, 343, 107 

Bell, E. F. 2003, ApJ, 586, 794 

Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 

Boselli, A., Gavazzi, G., Donas, J., & Scodeggio, M. 2001, AJ, 121, 753 

Brinchmann, J., & Ellis, R. S. 2000, ApJ, 536, L77 

Buat, V., & Xu, C. 1996, A&A, 306, 61 

Buat, V., & Burgarella, D. 1998, A&A, 334, 772 

Buat, V., Donas, J., Milliard, B., & Xu, C. 1999, A&A, 352, 371 

Buat, V., et al. 2005, ApJ, 619, L51 

Burgarella, D., Buat, V., & Iglesias-Paramo, J. 2005, MNRAS, 360, 1413 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Cardiel, N., Elbaz, D., Schiavon, R. P., Willmer, C. N. A., Koo, D. C, Phillips, A. C, & 
Gallego, J. 2003, ApJ, 584, 76 

Charmandaris, V., Le Floc'h, E., & Mirabel, I. F. 2004, ApJ, 600, L15 



-26- 

Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 

Dale, D. A., Helou, G., Contursi, A., Silbermann, N. A., & Kolhatkar, S. 2001, ApJ, 549, 
215 

Devriendt, J. E. G., & Guiderdoni, B. 2000, A&A, 363, 851 

Flores, H., et al. 1999, ApJ, 517, 148 

Garnett, D. R., & Shields, G. A. 1987, ApJ, 317, 82 

Gavazzi, G., Pierini, D., & Boselli, A. 1996, A&A, 312, 397 

Gordon, K. D., Calzetti, D., & Witt, A. N. 1997, ApJ, 487, 625 

Gordon, K. D., Clayton, G. C, Witt, A. N., & Misselt, K. A. 2000, ApJ, 533, 236 

Guiderdoni, B., & Rocca-Volmerange, B. 1987, A&A, 186, 1 

Heavens, A., Panter, B., Jimenez, R., & Dunlop, J. 2004, Nature, 428, 625 

Heckman, T. M., et al. 2005, ApJ, 619, L35 

Hirashita, H., Buat, V., & Inoue, A. K. 2003, A&A, 410, 83 

Hoessel, J. C, Abbott, M. J., Saha, A., Mossman, A. E., & Danielson, G. E. 1990, AJ, 100, 
1151 

Iglesias-Paramo, J., Buat, V., Donas, J., Boselli, A., & Milliard, B. 2004, A&A, 419, 109 

Jarrett, T. H., Chester, T., Cutri, R., Schneider, S., Skrutskie, M., & Huchra, J. P. 2000, 
AJ, 119, 2498 

Karachentsev, I. D., et al. 2002, A&A, 389, 812 

Kennicutt, R. C. 1998, ARA&A, 36, 189 

Kong, X., Chariot, S., Brinchmann, J., & Fall, S. M. 2004, MNRAS, 349, 769 
Leitherer, C, et al. 1999, ApJS, 123, 3 

Lilly, S. J., Le Fevre, O., Hammer, F., & Crampton, D. 1996, ApJ, 460, LI 
Lonsdale Persson, C. J., & Helou, G. 1987, ApJ, 314, 513 



-27- 

Lynden-Bell, D. 1971, MNRAS, 155, 95 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., Steidel, C. C, & Fruchter, A. 
1996, MNRAS, 283, 1388 ApJS, 107, 215 

Martin, D. C, et al. 2005a, ApJ, 619, LI 

Martin, D. C, et al. 2005b, ApJ, 619, L59 

Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 64 
Morrissey, P., et al. 2005, ApJ, 619, L7 

Moshir, M. et al. 1990, IRAS Faint Source Catalogue, version 2.0 

Panuzzo, P., Bressan, A., Granato, G. L., Silva, L., & Danese, L. 2003, A&A, 409, 99 

Peebles, P. J. E. 1993, Principles of Physical Cosmology (Princeton: Princeton University 
Press), pp. 119-121 

Rowan-Robinson, M., et al. 1997, MNRAS, 289, 490 

Salim, S., et al. 2005, ApJ, 619, L39 

Saunders, W., et al. 2000, MNRAS, 317, 55 

Sauvage, M., & Thuan, T. X. 1992, ApJ, 396, L69 

Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11, 1 

Schiminovich, D., et al. 2005, ApJ, 619, L47 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, M., & Pettini, M. 1999, ApJ, 
519, 1 

Sullivan, M., Mobasher, B., Chan, B., Cram, L., Ellis, R., Treyer, M., & Hopkins, A. 2001, 
ApJ, 558, 72 

Takeuchi, T. T., Yoshikawa, K., Ishii, T. T. 2000, ApJS, 129, 1 
Takeuchi, T. T., Yoshikawa, K., Ishii, T. T. 2003, ApJ, 587, L89 
Takeuchi, T. T., Buat, V., & Burgarella, D. 2005, A&A, 440, L17 



-28- 

Vijh, U. P., Witt, A. N., & Gordon, K. D. 2003, ApJ, 587, 533 
Wang, B., & Heckman, T. M. 1996, ApJ, 457, 645 
Witt, A. N., & Gordon, K. D. 2000, ApJ, 528, 799 
Wyder, T. K., et al. 2005, ApJ, 619, L15 
Zaritsky, D. 1993, PASP, 105, 1006 



This preprint was prepared with the AAS iATgX macros v5.2. 



-29- 



Table 1: Typical uncertainties of the NUV and FUV magnitudes as a function of the mag- 
nitude. 



AB Magnitude 


ff(NUV) 


cr(FUV) 


interval 


(mag) 


(mag) 


< 16 


0.01 


0.02 


16 ... 18 


0.03 


0.05 


18 ... 20 


0.09 


0.15 


20 ... 22 


0.26 


0.40 



Table 2: Basic properties of the sample galaxies: (1) Name; (2) Flag indicating membership to 
the NUV-selected subsample; (3) Flag indicating membership to the FIR-selected subsample; 
(4) R.A. (J2000); (5) Dec. (J2000); (6) Heliocentric velocity; (7) Distance derived from the 
velocity corrected for the Local Group infall onto Virgo and Ho = 70 km s -1 Mpc^ 1 ; (8) 
Morphological type from NED; (9) IRAS identificator: F for FSC origin; Q,0,R for PSCz 
origin; SCANPI for absence in both catalogs. 



Name 


UVsel 


FIRsel 


R.A. 


Dec. 


vel. 


Dist 


Type 


IRAS Id. 








(J2000.0) 


(km sec -1 ) 


(Mpc) 






MRK 544 


Y 


N 


4 48.70 


- 1 29 54.60 


7110 


101.39 


S? 


F00022-0146 


NGC 10 


Y 


Y 


8 34.56 


-33 51 27.25 


6811 


94.62 


Sbc 


Q00060-3408 


NGC 35 


Y 


Y 


11 10.46 


-12 1 14.74 


5964 


83.95 


Sb 


Q00086-1217 


NGC 47 


Y 


Y 


14 30.42 


- 7 10 6.28 


5700 


80.54 


Sbc 


Q00119-0726 


NGC 101 


Y 


N 


23 54.72 


-32 32 9.06 


3383 


45.92 


Sc 


F00214-3248 



Note: The distances to UGCA 438 and UGC 12613 were taken from Karachentsev et al. (2002) and Hoessel 
et al. (1990) respectively. 
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Table 3: Photometric properties of the sample galaxies: (1) Optical identifier; (2) NUV 
magnitude corrected for Galactic extinction; (3) FUV magnitude corrected for Galactic 
extinction; (4) Flux density at QOfim in Jy; (5) Flux density at 100/xm in Jy; (6) H magnitude 
from 2MASS Extended Source Caztalog. For galaxies with no detection at 2MASS we adopt 
the limiting value of H = 13.9 mag (3 mJy) as given by Jarrett et al. (2000); (7) NUV 
attenuation in mag, derived as indicated in Buat et al. (2005); (8) FUV attenuation derived 
as indicated in Buat et al. (2005). For some galaxies Eqs. 1 and 2 gave negative values of the 
NUV and FUV attenuations which is senseless. In fact, this is an artifact of the polynomial 
fitting used to derive a functional form for the attenuation in Buat et al. (2005). Throughout 
this paper they will be considered as zero. 



Name 


ABnuv 


ABpuv 




/ioo 


H 


Anuv Afuv 




(ma; 


g) 


(Jy) 


(mag) 


(mag) 


MRK 544 


15.75 


15.97 


0.49 


1.11 


12.50 


0.71 0.96 


NGC 10 


15.48 


15.98 


0.66 


2.97 


9.50 


1.59 2.10 


NGC 35 


15.55 


15.89 


1.31 


2.34 


11.91 


1.29 1.67 


NGC 47 


15.46 


15.97 


0.85 


2.57 


10.25 


1.02 1.49 


NGC 101 


14.78 


14.98 


0.55 


1.75 


10.74 


0.50 0.72 



Table 4: Star formation related properties of the sample galaxies: (1) Identifier of the galaxy; 
(2) SFR NUV from Eq. 3; (3) SFR FUV from Eq. 4; (4) SFR dust from Eq. 5; (5) SFR tot (NUV) 
from Eq. 6; (6) SFR tot (FUV) from Eq. 6 but using SFR° FUV instead of SFR° NUV ; (7) (SFR) 
averaged over the galaxy's lifetime, estimated as indicated in Appendix B. 



Name 


log SFR NUV 


log SFR FUV 


log SFR dust 


log SFR tot (NUV) 
(M yr- 1 ) 


log SFR tot (FUV) 


log (SFR) 


MRK 544 


0.83 


0.85 


0.60 


0.80 


0.75 


0.81 


NGC 10 


1.23 


1.24 


1.07 


1.09 


1.02 


1.95 


NGC 35 


0.98 


1.00 


0.78 


0.84 


0.79 


0.88 


NGC 47 


0.87 


0.86 


0.78 


0.85 


0.77 


1.51 


NGC 101 


0.45 


0.46 


0.13 


0.43 


0.38 


0.83 
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Fig. 1. — NUV (blue filled dots) and 60/im (red stars) LFs for the NUV and FIR-selected 
samples. The blue and red lines correspond to the NUV and 60/im LFs from Wyder et al. 
(2005) and Takeuchi et al. (2003) respectively. Error bars correspond to l-a uncertainty. 
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Fig. 2. — (a) Redshift distributions for the NUV (blue) and FIR (red) selected galaxies, (b) 
Theoretical redshift distributions for samples of galaxies with the same limiting fluxes as our 
samples, (c) Same as (b) but for a NUV-selected sample with AB^uy < 18 mag. Vertical 
dashed lines in the three panels correspond to the median values of each distribution. 
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Fig. 3. — L 60 vs. Lnuv (non corrected for attenuation) for the NUV (blue circles) and 
FIR (red stars) selected samples. Arrows indicate upper limits. Dashed curves represent 
the loci of points with L to t = 10 7 , 10 8 , 10 9 , 10 10 and 10 n L & respectively. The dotted line 
corresponds to L NUV = L 60 . The red (blue) dashed line corresponds to the maximum 
(minimum) attenuation below (above) which the FIR (NUV) selected sample is complete. 
The larger stars and circles correspond to the UVLGs of the FIR and NUV selected samples 
respectively. 
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Fig. 4. — log SFR NUV / SFR FUV vs. L tot . Symbols are as in Figure 3. 
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Fig. 5. — log S F R NUV / S F R dust vs. L tot . Symbols are as in Figure 3. The right Y axis 
indicates the NUV dust attenuation. 



-36- 



0.30 



-i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — r~ 



0.20 



> 



DC 
00 

\ 

;> 0.10 



Li- 
en 



o 



****** 



0.00 



-0.10 



_l I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I l_ 



-1 



log SFR (ot (NUV) 



Fig. 6.— log SFR tot (NUV)/SFR tot (FUV) vs. log SFR tot (NUV). Symbols are as in Fig- 
ure 3. 
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Fig. 7.- (a) log SFR NUV /SFR tot (NUV) vs. log SFR tot (NUV) (b) 

log SFR dnst /SFR t ot(NUV) vs. log SFR tot (NUV). Symbols are as in Figure 3. 
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Fig. 8. — Lm vs. Lm/Lpuv- Symbols are as in Figure 3. The dashed line corresponds to 
the relation L m /L FUV ~ y/L m /10 9 of Bell (2003). 



-39- 



40 



30 



o 
cn 

o 

- 20 
n 



10 



n i i i i i i i i n i i i r 



1 1 ' 



n i i i i i n 



-l l l I l J _l l I l l l l I l L 



l _ jfe&>>>_>_: 

_L 







-1.5 



-1.0 -0.5 



0.0 
log b 



0.5 



1.0 



1.5 



Fig. 9. — Histogram of the distribution of log b for the NUV (blue dashed line) and the FIR 
(red solid line) selected samples. The vertical dotted lines correspond to the median values 
of the distributions. Arrows correspond to nominal limiting values. 
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Fig. 10. — (a) SFRtot vs. b. (b) Anuv vs. b. Symbols are as in Figure 3. 
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Fig. 11. — H-band luminosity vs. SFR to t- Symbols are as in Figure 3. 
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Fig. 12. — if-band luminosity vs. A^uv- Symbols are as in Figure 3. 
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Fig. 13. — H-hand luminosity vs. b. Symbols are as in Figure 3. The dashed line corresponds 
to the calibration of Boselli et al. (2001) for their sample of late-type galaxies from nearby 
clusters with a normal Hi content. 



